Bound diquarks and their Bose-Einstein condensation 
in strongly coupled quark matter 
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We explore the formation of diquark molecules and their Bose-Einstein condensation (BEC) in 
the phase diagram of three-flavor quark matter at nonzero temperature, T, and quark chemical 
potential, /i. Using a quark model with a four-fermion interaction, we identify possible diquark 
excitations as poles of the microscopically computed diquark propagator. The quark masses are 
obtained by solving a dynamical equation for the chiral condensate and are found to determine the 
stability of the diquark excitations. The stability of diquark excitations is investigated in the T — fi 
plane for different values of the diquark coupling strength. We find that bound diquark molecules 
appear at small quark chemical potentials at intermediate coupling and that BEC of non-strange 
diquark molecules occurs if the attractive interaction between quarks is sufficiently strong. 



I. INTRODUCTION 

The one-gluon exchange interaction between quarks is 
attractive in the color-antitriplet channel and leads to 
color superconductivity in cold and dense quark matter 
Due to asymptotic freedom of QCD, the strength of 
the gluon-exchange interaction around the Fermi surface 
changes with the quark number density. At asymptot- 
ically high densities, the interaction is sufhciently weak 
to apply perturbation theory and the mean-field approx- 
imation, and color superconductivity is very similar to 
standard BCS superconductivity. The gap parameter is 
A ~ /xexp(— 1/(7), where g is the QCD couphng con- 
stant. Quark Cooper pairs are large, with a correlation 
length ^ ~ 1/A which is parametrically larger than the 
interparticle distance d ~ On the other hand, at 

intermediate densities which may be realized in the cores 
of compact stars and/or in the intermediate stages of 
heavy- ion collisions, the quark-quark interaction is rel- 
atively strong and the properties of the quark Cooper 
pairs will be modified. In particular, we expect their size 
to become of the same order as the interparticle distance, 
^ ~ 1/A ~ d. Moreover, it was argued that, due to the 
strong coupling, the fluctuations of the diquark-pair field 
become large around the critical temperature ,.2:j and 
that they give rise to a pseudogap region in the normal 
phase above Tc [1, 3 • 

If the quark-quark interaction becomes strong enough 
before the quarks are confined at lower density, quark 
Cooper pairs become small enough to be considered 
as diquark molecules which are tightly bound states of 
two quarks. The color-superconducting ground state 
of quark matter at low temperature turns into a Bose- 



Einstein condensed phase of diquark molecules 0, 0]. 
Most likely, the transition happens continuously with the 
change of the density, just like the BCS-BEC crossover 
in conventional condensed matter systems 0, § , see also 
Refs. [i,|ll[illlj,ll3|,lll[ll. An interesting property of 
diquark molecules is that they can exist even above the 
critical temperature Tc of BEC, up to the dissociation 
temperature Tdiss > [§]. If such modes exist in quark 
matter, they will affect its properties and experimental 
observables. 

In this work, we explore the appearance of diquark 
molecules and their BEC in the phase diagram of quark 
matter using a low-energy effective model. This model 
features an attractive quark-quark interaction with a con- 
stant coupling strength Gd that is regarded as a free pa- 
rameter of the model. We show that diquark molecules 
appear at low density at intermediate values of Gd. It is 
also shown that BEC of diquark molecules can occur for 
large values of Gd. 

In the normal phase above Tc, the strongest decay 
mode of diquarks is that into two quarks. Since the exci- 
tation energy of a quark at rest is M — /i, with M being 
the mass of the quark, the threshold energy for this de- 
cay process is (jJthr = 2(M — p.), where M and ft are the 
average mass and chemical potential of the quarks in the 
diquark. Recalling that the energy of the diquark exci- 
tations should be positive to ensure the stability of the 
system, one finds that the necessary condition for the 
existence of stable diquarks is Wthr > 0, or 



M > p,. 



(1) 
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As we will see later, a.t T = Tc Eq. ([T]) is also a suffi- 
cient condition [1, Q . From Eq. ([J) , we conclude that 
the stability of diquark excitations is determined by the 
quark masses. Quark masses are dynamically generated 
by chiral symmetry breaking and change as functions of 
T and 11. In this work, we solve the gap equations for 
the chiral condensates and incorporate this effect in our 
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calculation. Equation ^ also indicates that diquarks 
composed of heavier quarks tend to be more stable than 
those composed of light quarks, provided they exist at 
all. One thus expects that diquarks including a strange 
quark are more stable than those composed of up and 
down quarks. 

Usually, BEG is discussed using the canonical ensem- 
ble, i.e., the particle number density is fixed as an ex- 
ternal parameter. In this work, however, we employ the 
grand canonical ensemble and draw the phase diagram in 
the T — fi plane, as is usually done in the literature when 
exploring the QCD phase diagram T|. In order to decide 
whether BEG of diquarks occurs in this ensemble, we re- 
gard the region of the superconducting phase satisfying 
Eq. Il]) as Bose-Einstein condensed phase [sj]. 

In this exploratory study, we employ a common chem- 
ical potential /i for all flavors and colors. For quark mat- 
ter in compact stars, this is probably not a very good 
assumption, as the chemical potentials should be deter- 
mined to satisfy the neutrality and beta-equilibrium con- 
ditions. It is known that a rich phase structure can ap- 
pear under these conditions However, as we shall 
see in the following, diquark excitations play an impor- 
tant role even at high temperatures and small chemical 
potentials in the range relevant for heavy-ion collisions. 
In this case, our assumption of equal chemical potentials 
for all quark flavors and colors is applicable to very good 
approximation. 

This statement warrants a few remarks. Thermal 
model fits of hadron yields at chemical freeze-out show 
that, because of strangeness and isospin conservation, 
neither the strangeness chemical potential /zg nor the 
isospin chemical potential /x/ are zero. The strangeness 
chemical potential is nonzero because of associated pro- 
duction channels in hadronic matter at nonzero baryon 
chemical potential. For quark matter, however, fj,s 
should be strictly zero if the system has zero strangeness, 
and we may neglect fis- The isospin chemical potential 
is nonzero because of the initial isospin asymmetry of the 
colliding nuclei. However, at all collision energies it has 
been demonstrated that /x^ = 3/i Thus, to 

leading order we may set /i/ = 0. 



II. FORMALISM 

In order to study the phase diagram of quark matter, 
we employ a three-flavor quark model with four-fermion 
interactions. The Lagrangian is given by 

8 

a=0 



7,c 



(2) 



where the quark field tpa has color (a = r, g, b) and flavor 
(a = M, d, s) indices. The matrix of current quark masses 



is given by m = diag^(m„, md, m^); Aq = ^/2/3t and 
Aq, a = 1,...,8, are the Gell-Mann matrices in flavor 
space. The charge conjugate spinors are ipc = C?/'^ and 
ipQ = ip'^C, where C = 17^7'^ is the charge conjuga- 
tion matrix. In the following, we only consider diquark 
condensates and diquark excitations in the color anti- 
triplet channel. For the numerical calculations, we em- 
ploy a three-dimensional momentum cutoff A. We treat 
the diquark coupling constant G^) as a free parameter. 
For the other parameters, we use the values of Ref. p^ . 
TO„ = mrf = 5 MeV, m, = 120 MeV, Gs = 6.41 GeV"^ 
and A = 600 MeV. 

We evaluate the thermodynamic potential in the mean- 
field approximation; 



4Gd 



c=l 



T 
2 



/cPp 
^^-p-Tro.fx In [S-^{iujn,p)] 



(3) 



where the trace is taken over Dirac, flavor, and color 
indices, cun — (2n + l)n/T is the Matsubara frequency 
for fermions, and 



Ma = ma - AGsi^pall^o 

A,^2Gd{i^cPc^), 



(4) 
(5) 



are the constituent quark masses and the gap parameters 
for color superconductivity, respectively, with {Pc)'^p — 
ij5e°'^'^eabc- The 72 x 72 Nambu-Gor'kov propagator is 
defined by 



S ^{iuJn,p) 



(6) 



The 



with = iujnjo -Pl and M = diag ^ {M^ , Md , 
quark chemical potential /i has a common value for all 
flavors and colors, as mentioned above. 

The physical values of the variational parameters Ac 
and Ma satisfy the stationary conditions (the gap equa- 
tions) 



dn on 

and 



dA 



dMa 



0. 



(7) 



As we will see later, the color-superconducting phase 
transitions at nonzero temperature are of second order. 
Thus, the critical temperatures are determined by solving 
the following equation: 



1 



Ac 9Ac 



0. 



(8) 



We shall see later that this equation determines the poles 
of the diquark propagator at vanishing energy and mo- 
mentum. 

Since up and down flavors are degenerate in our model, 
we always have M„ = Md and Ai = A2. In the following. 
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we refer to the phase with A3 ^ and A1.2 = as the 
2SC phase; the phase with A3 ^ and Ai^2 7^ is 
the CFL phase Unpaired quark matter has Ai = 
A2 = A3 = 0. Because of the exphcit chiral symmetry 
breaking by the nonzero current quark masses, the chiral 
condensates (i^itpi) never vanish. 

At nonzero temperature, the order parameter fields, 



Ac{x,t) and Ma{x,t) have fluctuations around the val- 
ues determined by the mean-field approximation. In 
the following, we consider the amplitude fluctuations of 
Ac{x,t) in unpaired quark matter. The propagation of 
the fluctuations is characterized by the retarded propa- 
gator 



D^ix, t- x', t') = -ie{t - t'){[{l;{x, t)P,i,cix, t), i^cix', t')P,i^{x', t')]) = 



(Pkdu 



i') zp-(x— x') 



(9) 



where c = 1,2, and 3 correspond to the down-strange, up-strange, and up-down diquark fields, respectively. In the 
random-phase approximation, the diquark propagators are given by 



1 



2 1-|-Gz5Q«(p,w 

where Q^{p,Lu) is the one-loop quark-quark polarization function. In the imaginary-time formalism, it is given by 



(10) 



/d q 
^^--^TrD,f,c[^'c5'o(p - q,ii^n - ii^m)PcC {q,iuJm)C], 



(11) 



where Vn = 2Trn/T denotes the Matsubara frequency for bosons, and SQ{p,iujn) = + M7o — M] ^. Taking the 
analytic continuation, (5^(p, w) — Qc{p, 'ii^n)\w„->^+iv' obtain 



^ 1 V / c i 1 



s,t=± 



{sEp + tE^Y - IpP - {6M,f fjtE^ - /i) - f{-sE0 + ^l) 
stEpE^ uj + 2fi — sEp — tE^ + irj ' 



(12) 



where Ep = ^\q-p\^ + M^, E-y = ^\q\^+M^, SMc = \Mp - Af^|, and j{E) = [cxp(i;/r) + 1]-^ is the Fermi 

distribution function. The imaginary part of Q^{p,uj) denotes the difference of decay and production rates of the 
diquark field. At p = 0, it is given by 



><{-[(!- //3^)(1 - fy) - f^f^] + - - ^7) + [(1 - fp)i^ - f^) - fpf^] + 2^i + Ep + Ey) 
- fp{l - f^) - (1 - fp )f^] + 2^i + Ep^ E,) - [/+(1 - /-) - (1 - /+)/-] 5{uj + 2^i-Ep + E.,)] 



(13) 



where /± = {exp[(£;„ T m)/?"] + 1}"^- The first (second) 
term in the bracket in Eq. psp includes the decay pro- 
cesses of the diquark into two quarks (anti-quarks) and 
takes nonzero values at w > 2Mc — 2fj,{u!< —2Mc — 2fi), 
with Mc = {Mp + Mj)/2. The third and fourth terms 
represent Landau damping of the diquark. These terms 
become nonzero at —SMc — 2)1 < lo < SMc — 2fi. 

The poles of the diquark propagator are deter- 
mined by solving D^{p,uj) — 0, or equivalently 

l + Gz5Qf(p,c^) = 0, (14) 



in the complex-energy plane. Applying lu = |p| = to 
this equation, one can easily show that Eq. (|14|) is equiva- 
lent to the critical condition Eq. ([8]). Therefore, -D^(p, w) 
has a pole at the origin at T = Tc of the second order 
transition. This property is known as the Thouless cri- 
terion in condensed matter physics Above the 
pole moves continuously from the origin to the fourth 
quadrant. This mode is called the soft mode. If Mc < /i 
at T = Tc, w = is in the continuum of the decay pro- 
cess into two quarks and the imaginary part at the pole 
starts growing just above Tc 0]. When Mc > /i, on the 
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other hand, the soft mode is stable against spontaneous 
breaking into a pair of quarks and the pole moves on the 
real axis in the vicinity of Tc- This mode is nothingbut a 
bound state of two quarks: the diquark molecule 8]. As 
T increases, the pole eventually arrives at the threshold 
of the decay process into two quarks Wthr = 2(Mc — at 
the dissociation temperature Tdiss, and the soft mode is 
no longer a bound state at T > T^j^^. Since the pole is at 
wthr at Tj^iss, the dissociation temperature is determined 
by solving 

1 + GdQc{P = 0,Wthr) I T=T,-,^^ = 0. (15) 

Although the diquark modes can acquire decay rates 
due to Landau damping, i.e., the third and fourth term 
in Eq. ((T3|), our numerical results show that the soft 
modes never appear in the range of energies where Lan- 
dau damping is nonzero. Therefore, these processes do 
not contribute to the decay rate of the soft modes in the 
parameter range employed in the present study. There 
can appear another pole of D^{p,uj) instead of the soft 
mode at the energy —2Mc — 2fjL < oj < ~5Mc — 2/x. This 
mode does not have a decay rate and should be identified 
as a bound anti-diquark Q. For lower /z, thermal excita- 
tions of bound anti-diquarks play an important role. 

If bound diquarks are formed at T = Tc, it is natural 
to identify the color-superconducting phase below Tc as 
a Bose-Einstein condensed phase of diquark molecules 
Q. In the following, therefore, we regard the color- 
superconducting phase satisfying Mc < /i as a Bose- 
Einstein condensate. Notice, however, that this is just 
a rough guide to separate the BEC and BCS regions; 
these two limits are connected continuously and there is 
no sharp phase boundary between them [Sj. 

III. NUMERICAL RESULTS 

In this section, we show the phase diagram in the T-/i 
plane for several values of the diquark coupling Gd- In 
Fig. [11 we first discuss the case Gd/Gs = 0.75, which 
is a value commonly used in the literature [J, hj. The 
bold and thin solid lines represent first- and second-order 
phase transitions, respectively. One sees that there ap- 
pear two types of color-superconducting phases, the 2SC 
and CFL phase, at high /i and low T. At T = 0, these 
phases are separated by a first-order phase transition; 
the first-order transition terminates at a critical point 
for some nonzero value of temperature. The dissocia- 
tion temperatures of diquark molecules T^^^^s ^^^^ shown 
by the dashed lines. We see that there exists a region at 
small chemical potential where stable diquark molecules 
are formed. 

In order to see whether diquark molecules undergo 
BEC, we plot the conditions /i ~ Mc by the dash- 
dotted hues in Fig. [U the regions to the left of these 
lines satisfy ^ < Mc- We see that these lines termi- 
nate at the first-order transition and we do not have a 
color-superconducting phase satisfying fi < Mc- In other 
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FIG. 1: The phase diagram in the T-p plane for Gd/Gs = 
0.75. The bold and thin solid lines represent first- and second- 
order phase transitions. The dashed lines denote the dissocia- 
tion temperature of diquark molecules for up-down diquarks, 
T^igj,, and for up-strange and down-strange diquarks, T^(^^- 
The conditions n = Mc are shown by the dash-dotted fines. 



words, BEC does not appear for this value of Gd/Gs- 
The behavior of Mc and Ac as functions of at T = 
are shown in the upper panel of Fig.[2l One observes that 
M^ = Mu^ has a discontinuity at ~ 343 MeV corre- 
sponding to a first-order transition, and M3 is larger than 
^ to the left of the discontinuity. The diquark conden- 
sate A3 assumes nonzero values only for /i > M3. This 
is a typical property at weak coupling; when /i > M3, 
the Fermi surfaces of up and down quarks exist and the 
Cooper instability leads to a diquark condensate, while if 
not, the ground state is nothing other than the vacuum. 

It is worth mentioning that bound diquark molecules 
appear in the phase diagram even though BEC does not 
exist in the phase diagram. The diquark coupling used 
in Fig. [1] is strong enough to form bound diquarks, but 
it is still too weak to lead to their BEC. 

Next, we show the phase diagrams with much stronger 
diquark couplings. In Fig. [3l the phase diagrams with 
Gd/Gs = 1-1 and 1.2 are shown. We see that, as Gd 
becomes larger, the regions of the 2SC and CFL phases 
expand toward lower /i and higher T. For Gd/Gs = 1-1, 
there appears BEC of up-down diquarks in the region 
of the 2SC phase satisfying fi < M3, shown by the 
shaded area in Fig. [3] The BEC region becomes wider 
for Gd/Gs = 1.2. One also observes that the dissocia- 
tion temperatures T^j^g become higher as Gd increases. 
In the phase diagrams in Fig. [31 T^^^^ at fi — are com- 
parable or much higher than the critical temperature of 
the QCD phase transition, which is predicted to be in the 
range Tc ^ 150 — 190 MeV in lattice QCD simulations 
[2O I. This result shows that diquark molecules can ex- 
ist even in the quark-gluon plasma phase if the diquark 
coupling is strong enough. 
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FIG. 2: Order parameters Ma and Ac at T = as functions 
of /i for various values of the diquark coupling Gd/Gs ~ 
0.75, 1.1 and 1.2. The chemical potential is shown by the 
dashed line. 



In order to see the diquark coupling dependence of 
the dissociation temperatures, we show T^-^gg at = 
as functions of Gd/Gs in Fig. [4l At weak coupling, 
^diss — ^ bound diquarks do not exist. As Gd/Gs 
becomes larger, T^-^^^ eventually become nonzero and in- 
crease rapidly. The dissociation temperatures for di- 
quarks including the strange quark, T^l^^, are always 
higher than that for the up-down diquark. This fea- 
ture comes from the difference of the threshold energy 
2(Mc-m)- 

The other interesting feature shown in Figs. [1] and [3] is 
the behavior of the line of first-order phase transitions. 
The first-order phase transition at lower density is shorter 
for Gd/Gs = 1.1 than that for Gd/Gs = 0.75, and dis- 
appears at Gd/Gs = 1.2. To understand this behavior, 
we display the order parameters for Gd/Gs = 1.1 and 
1.2 in the middle and lower panels of Fig. [2l respectively. 
The figure shows that, as Gd becomes larger, Ac increase 
while the quark masses Mc become smaller, and the dis- 
continuity of M3 disappears at Gd/Gs — 1.2. The de- 
crease of Mc can be understood as the interplay between 
the chiral and diquark condensates [2l|: the energy gain 
due to diquark condensation is proportional to the sur- 
face area of the Fermi sphere. Since the radius of the 
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FIG. 3: The phase diagram in the T-fi plane for relatively 
strong diquark couplings Gd/Gs = 1.1 and 1.2. BEG of up- 
down diquarks occurs in the shaded area. 
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FIG. 4; Dissociation temperatures of up-down diquarks T^ig, 
and up-strange and down-strange diquarks T^^^^- 
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FIG. 5: The phase diagram in the T-jj, plane for extremely 
strong diquark coupling Gd/Gs = 1-5. 



Fermi sphere is pp = a/zz^ — M^, the condensation en- 
ergy increases when M decreases. Since the masses of 
quarks are suppressed as Gjj becomes larger, the lines 
for /i = Ale and the region of BEC move toward lower /i. 

Finally, let us consider the phase diagram for extremely 
large diquark coupling. In Fig. [51 we show the phase 
diagram for Gd/Gs — 1-5. The 2SC and CFL phases 
become much wider toward lower /x and higher T, and 
T^igg becomes higher. The region of BEC of up-down 
diquarks also becomes wide. We do not obtain a region 
in the CFL phase satisfying /x > Mi_2, i.e., BEC of up- 
strange and down-strange diquarks does not occur for 
Gd/Gs < 1-5. If the diquark couphng is raised further. 



the vacuum, i.e., T = = 0, eventually becomes a Bose- 
Einstein condensate of diquark molecules, which is clearly 
unphysical. 



IV. SUMMARY AND DISCUSSIONS 

In this paper, we explored the phase diagram of three- 
flavor quark matter focusing on the appearance of di- 
quark molecules and their Bose-Einstein condensation 
under variation of the diquark coupling constant Gd- 
We found that diquark molecules can appear at small 
/i and (probably realistically large) intermediate values 
of the diquark coupling, while BEC of up-down diquarks 
is realized for (probably unphysically) large values of the 
diquark coupling. The dissociation temperatures of di- 
quarks become higher as Gd increases. At strong cou- 
pling, the dissociation temperatures could be higher than 
the critical temperature of the deconfinement transition. 

In this work, we employed the random-phase approx- 
imation for the calculation of the diquark propagator. 
In this approximation, the propagators in the polar- 
ization function Eq. (jlip are those for non-interacting 
quarks, i.e., the effect of diquark excitations is not self- 
consistently incorporated. An extension of the present 
work to include this effect is an interesting subject for 
further study, because we expect that the formation of 
diquark molecules would modify the result especially at 
high temperatures. Similarly, the mean-field approxima- 
tion used to draw the phase diagram is no longer appli- 
cable due to the existence of well-developed soft modes 
in a strongly coupled system. The incorporation of these 
effects has already been partially made in Ref. 15]. 

M.K. thanks H. Abuki and T. Kunihiro for discussions. 
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